Atomic Collapse and Quasi-Rydberg States in Graphene 
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Charge impurities in graphene can host an infinite family of Rydberg-like resonance states of 
massless Dirac particles. These states, appearing for supercritical charge, are described by Bohr- 
Sommerfeld quantization of collapsing classical trajectories that descend on point charge, in analogy 
to Rydberg states relation with planetary orbits. We argue that divalent and trivalent charge 
impurities in graphene is an ideal system for realization of this atomic collapse regime. Strong 
coupling of these states to the Dirac continuum via Klein tunneling leads to striking resonance 
effects with direct signatures in transport, local properties, and enhancement of the Kondo effect. 



The discovery of massless Dirac excitations in 
graphene [H triggered new interest in solid-state realiza- 
tion of quantum electrodynamics (QED) 0, S, 0] • Trans- 
port phenomena in this system p, 0] can be used to 
probe classic concepts of QED, such as chiral dynam- 
ics Q, flavor degrees of freedom Q and particle/hole co- 
existence Q . Here we demonstrate that graphene opens a 
way to investigate in the laboratory a fundamental quan- 
tum relativistic phenomenon, that is, atomic collapse in 
a strong Coulomb electric field [7J, LLO( , long sought for but 
still inaccessible in high-energy experiments 11 1 . 

Bohr's theory of an atom has explained that, while 
an electron is irresistibly pulled to the nucleus by the 
Coulomb force, it is prevented from falling on it by the 
quantum mechanical zero-point motion. This balance, 
however, becomes more delicate in the relativistic the- 
ory. The effects undermining the stability of matter arise 
already in classical dynamics, where electron trajectory 
can spiral around the nucleus and eventually fall down 
on it {l2j| (see FigQJt,b), provided that electron angular 
momentum is small enough: M < M c = Ze 2 /c, where 
Z is nuclear charge. Quantum mechanics partially saves 
matter from collapse by imposing the angular momen- 
tum quantization M = nh, which makes the relativistic 
fall-down possible only for heavy Z > hc/e 2 w 137. 

Early work on the Dirac-Kepler problem has revealed 
bizarre properties of atoms with nuclear charge in excess 
of Z = 137, posing as a fundamental bound on the peri- 
odic table of elements extent at large Z. The breakdown 
at Z > 137 of the \ow-Z solution of the Dirac equation 
requires accounting for a finite nuclear radius [231 ] . The 
resulting electron states dive into the hole continuum at 
Z > 170 and decay by positron emission 
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phenomena, never observed in the laboratory due to the 
difficulty of producing heavy nuclei, should be more read- 
ily accessible in graphene owing to its large "fine struc- 
ture constant," a = e 2 /hvF ~ 2.5, where vp w 10 6 m/s 
is the velocity of Dirac excitations. 
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FIG. 1: Classical and quantum picture of atomic collapse due 
to electron with angular momentum M < M c falling on the 
nucleus. Trajectories obtained from {1} for (A) positive and 
(B) negative energy e of a massless Dirac particle are shown. 
(B,C) At e < there are collapsing particle trajectories and 
non-collapsing hole trajectories, separated by a classically for- 
bidden region, the annulus n < r < ri (ri,a = r, =p Mvf/e, 
r» = Ze 2 /|e|). Tunnel coupling to the continuum at r > T2 
defines a family of quasistationary states with complex en- 
ergy spectrum (D), appearing abruptly when the potential 
strength exceeds the stability threshold \/3\ = |. 



and is thus regarded as one of the main factors limiting 
carrier mobility. Recent inv estig ations of screening of im- 



Charge impurities are an essential ingredient of our 
current understanding of transport in grap hene. Scatter- 
ing on charge impurities explains the linear 
density dependence of conductivity in this material [l6| , 



purity potential 17, 3, 3 2oT 21| have only reinforced 
these conclusions, making graphene an ideal test system 
for the theory of Coulomb scattering [3, 0] of massless 
Dirac particles. 

In this work we show that, although massless particles 
cannot form bound states, an infinite family of quasi- 
bound states appears abruptly when the Coulomb po- 
tential strength exceeds a certain critical value j3 = \. 
These states are obtained from Bohr-Sommerfeld quan- 
tization of collapsing orbits which descend on the point 
charge, similar to how the hydrogenic Rydberg states are 
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found from circular orbits. The energies of these states 
converge on zero, e„ — > at large n, whereas their radii 
diverge, similar to the Rydberg states. These results are 
corroborated by an exact solution of the 2D Dirac-Kepler 
problem. In graphene, the effective Coulomb potential 
strength is given by (3 — Ze 2 / nhvp with intrinsic di- 
electric constant k w 5, and therefore the critical value 
(3 = i can be reached already for the impurity charge 
Z > 1. This is a lot more convenient from the experi- 
mental point of view than Z > 170 in heavy atoms. 

Coupling of quasi-Rydberg states to the Dirac contin- 
uum, mediated by Klein tunneling 0] , leads to strong res- 
onances in the scattering cross-section, manifest in trans- 
port, and to striking effects in local properties that can 
be probed by deliberately introducing impurities with 
Z > 1 in graphene. Univalent charge impurities, such 
as K, Na, or NH3, all commonly used in graphene, are 
on the border of the supercritical regime. To investigate 
this regime experimentally, one can use divalent or triva- 
lent dopants such as alkaline-earth or rare-earth metals. 
They are frequently used to prepare intercalated graphite 
compounds [26j|, e.g., Ca and Yb [27] (Z = 2), La [28( and 
Gd [29} (Z = 3). Recently, spectroscopic experiments on 
graphene doped by Ca have been reported [30] . 

Recent literature [13, EE E H E| investigated the 
problem of screening of charge impurities, which depends 
on the polarization of the Dirac vacuum [18j . This effect 
is mostly inconsequential in atomic physics due to its 
short spatial scale set by the Compton wavelength A = 
h/mc « 2.4 x 10~ 3 nm. In the massless case of graphene, 
however, it leads to long-range polarization, appearing 
above the critical value (3 = | of the impurity charge [l8| . 
These studies indicate that (3 = \ separates two very 
different regimes of screening, essentially perturbative at 
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20j, and nonlinear at (3 > 
To explain why the quasi-bound states appear at large 
f3, we consider fermions with energy e < in the potential 
V(r) = —Ze 2 jr. Since the kinetic energy K = e — V(r) 
vanishes at r» = Ze 2 /\e\, the polarity of carriers changes 
sign inside the disk r < r* (see FigQJ). If r* exceeds 
particle wavelength A = Hvf/\s\, which happens for [3 = 
Ze 2 /Hvf ^ 1, quantum states can be trapped at r < 
r*. These states will have finite lifetime due to Klein 
tunneling through the barrier at r ss r* . Crucially, since 
the ratio r*/A is independent of e, this reasoning predicts 
infinitely many quasi-bound states (see FigfTJi). 

These states can be constructed quasiclassically, from 
relativistic dynamics described by the Hamiltonian H = 
vf\p\ + V(r), where V(r) = —Ze 2 jr. The collapsing tra- 
jectories with angular momenta M < M c — Ze 2 /vf are 
separated from non-falling trajectories by a centrifugal 
barrier. This is manifest in the radial dynamics 
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sically forbidden region, the annulus r± < r < r2, 
r i,2 = (Ze 2 =F Mvp)/s, where ^ is negative. The quasi- 
bound states trapped by this barrier can be found from 
the Bohr-Sommerfeld quantization condition f 1 p r dr = 
irhn, where ro is the lattice cutoff (cf. Refs.0, [25J]). 
Evaluating the integral with logarithmic accuracy, we ob- 

Ze 2 _ „,u™„ „, — ( h,t2 _ ^2\ 1 / 2 
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irhn, where 7 = (M 2 



gives the quasi-Rydberg states 
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The energies @ are equally spaced on the log scale with 
the separation diverging as I/7 at the threshold M c w M. 

To find the transparency of the barrier, we integrate 
Inip r to obtain the tunneling action 



S : 
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Vf 



tt(M c - 7 ). (3) 



Taken near the threshold 7 ~ 0, the transparency e~ 2S l h 
gives the width r„ ~ |e n | exp(— 2irZe 2 /Hvf)- Notably, 
since S has no energy dependence, all the states © fea- 
ture the same width-to-energy ratio. 

It is instructive to compare these results to the exact 
solution of the Coulomb scattering problem. For that, 
we consider the Dirac equation for a massless electron in 
a potential V(r) — where (3 = —Ze 2 /hvF- (It will 
be convenient to include the minus sign in (3 to explic- 
itly account for attraction.) Performing standard angu- 
lar decomposition and solving the radial equation sepa- 
rately in each angular momentum channel [18], one finds 
the scattering phases S m (k) that behave differently in the 
three regimes: (i) /3 > (3 m , (ii) — f3 m < (3 < (3 m and (iii) 
(3 < -f3 m ((3m = \m + § |). As illustrated in Fig[2K, 
5 m (k) are energy-independent in case (ii) and have a log- 
arithmic dependence 5 m (k) ~ — 7mfcro in case (i), where 
7 = \/ ft 1 — Pm- I n the case (iii) the dependence is de- 
scribed by kinks of height tt equally spaced on a log scale: 



e 2zS m (k) = e 7ri/3„ 
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(see [18|, Eq.(20)) where 



e^T r(l+2ry) r(l-i 7 +z/3) 
Tj T(l-2i7) r(l+i7+i/3) ' 



(4) 



and 



where p r is the radial momentum. This defines a clas- 



X (k) = 2 7 ln2fcr + 2arctani±2, V=\/f^- (5) 

The average winding rate of the phase in case (iii), 
5 m {k) ~ 7 In Zero, is the same as in case (i) up to a sign 
(see Fig 12]). The kinks signal the appearance of quasi- 
bound states at negative energies. 

To find the quasi-bound states, we seek a scattering 
state with complex energy in which there is no incoming 
wave. This implies vanishing of the numerator (denom- 
inator) of ([3]) at e < (e > 0). In the e < case we 
obtain an equation for k: e lx ^ = —z. The right-hand 
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FIG. 2: (A) Scattering phase So at negative energy e — 
—VFk < 0. The kinks correspond to the quasi-bound states 
trapped by the impurity potential of supercritical strength 
ft < —1/2, as illustrated in Fig[T] (B) Transport cross- 
section (|10p vs. potential strength. Fano resonances cor- 
responding to quasi-bound states occur at ft < -i. The 
oscillatory behavior at ft > \ results from the energy de- 
pendence S m (k) ~ — 7 m lnfcro. The cross-section asymmetry 
upon ft — » — /3, with the values at ft < typically lower than 
at ft > 0, reflects that the Klein barrier prevents particles 
from reaching the region of strong scattering r ~ 0. 

side of this equation in general has a non-unit modulus, 
which makes it impossible to satisfy it by a real k. Com- 
plex solutions of = — zln(— z) — 2nin resemble those 
obtained quasiclassically, Eq.©. 

For a more direct comparison, let us consider (3 near 
the threshold (3 = (3 m . Expanding in small 7, we find 
solutions similar to our quasiclassical result |J5J| , 

C I 7T \ 

k n = exp n — iX ) , n > 0. (6) 

2^o V 7 / 

where A = 1 _ e _ 2j -g and the prefactor c is of order one. 
(We suppress n < 0, since Eq.© holds only for kr <1.) 

The interpretation of the solutions (O depends on the 
sign of p. For near-critical negative values « — \ we 
have argfc„ « 0.0457T, such that k n have small imaginary 
parts, defining sharp resonances of width 

ir= -Ime= ^\e\ «0.14|e|. (7) 

In contrast, argfc„ w — 1.0457T for positive (3 ps \, i.e. 
complex fc„'s are rotated by more than 180° away from 
positive semi-axis. Thus there are no long-lived states 
with positive k (s < 0). Instead, since argfc„ « 7r, in 
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FIG. 3: Peak structure in the ohmic conductivity (|8]). (I10[) 
for overcritical ft occurring when the Fermi energy is aligned 
with resonances. The values of ft are indicated near each 
trace. (Parameters used: ni mp = 3 ■ 10 11 cm" 2 , ro = 0.25 nm) 
Inset: Fano resonance structure in the cross-section (|10[) at 
negative energies with the overall 1/k dependence factored 
out. 



this case all fc„'s are found near the negative real semi- 
axis (e > 0). This is in agreement with particle/hole 
symmetry. 

Before discussing manifestations in graphene, where 
at finite carrier density the 1/r potential is screened, 
we note that the essential physics will be unaffected 
by screening as long as the quasi-bound states persist. 
At finite density, the RPA screening length is compara- 
ble to the Fermi wavelength Af = Tivf/£f [3, O EH, 
whereas our quasiclassical estimate of the state radius 
gives r\ « (M — M c )vf/\£\- The latter is much smaller 
than Xp near (3 = 4, which means that RPA screen- 
ing is non-detrimental for these states. Similarly, esti- 
mates for nonlinear screening [T3| indicate that its effect 
is inessential at weak coupling, leaving enough room for 
quasi-bound states. 

Resonance scattering on the quasi-bound states will 
manifest itself in the dependence of transport properties 
on the carrier density. Here we analyze electrical conduc- 
tivity described by the Drude-like model (see Ref.fpi]): 

e 2 

a = -j^Zzft, t 1 = v F n imp a tI , (8) 

where Ef is the Fermi energy, rii mp is the concentration 
of charge impurities and a tT is the transport scattering 
cross-section for one impurity. We use the 2D scattering 
amplitude partial wave decomposition 

r, ■ 00 

/fa) - -= ]T (e 2 ^ - 1) cos(m + ±>, (9) 

v m— 
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(see Ref.Q, Eq.(47)) to evaluate transport cross-section 



f 4 00 

<r tr = / d(f(l - cos ip)\f{if)\ 2 = - sin 2 ft 

^ m=0 



(10) 



0m = <5 m — d m +i, with the phases <5 m given by ([4]) for over- 
critical channels (see Ref. 22j for subcritical channels). 

For subcritical potential strength the phases are 
energy-independent and thus <r tr scales as l/|e|, giving 
conductivity ([5]) linear in the carrier density 13, 1J, 15| ]. 
For \j3\ > \ the contribution of the subcritical channels 
still scales as l/|e|, while the overcritical channels, be- 
cause of energy-dependent 6 m (k), give an oscillatory con- 
tribution (c/. Figl33). These oscillations, shown in Fig[3] 
inset, have a characteristic form of Fano resonances cen- 
tered at e„. In this regime the conductivity ([8]) exhibits 
peaks at the densities for which the Fermi energy £p 
alignes with e„. As evident from Fig[31 the peak position 
is highly sensitive to the potential strength /3, changing 
by an order of magnitude when [3 varies from —1.0 to 
— 1.3, which is a combined effect of Ef quadratic depen- 
dence on density and of the exponential dependence in 

Another striking feature in the conductivity vs. den- 
sity plots in FigJ3] is the n — ► — 11 asymmetry, which 
results from the scattering cross-section being typically 
lower at e < than at e > 0. Such asymmetry, noted al- 
ready in the subcritical regime [22| , becomes more promi- 
nent in the supercritical regime because of the Klein 
barrier preventing particles with negative energies from 
reaching the strong scattering region r ~ 0. 

The signatures of quasi-bound states, similar to those 
in conductivity, will be featured by other transport coef- 
ficients. In particular, they will be strong in the thermo- 
electric response because it is proportional to the energy 
derivative of Otr- Yet the most direct way to observe 
these states is via the local density of states (LDOS) 



mke,r)\\ k e 
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where ip is the two-component Dirac wave function 
(5) [1 81 ] . This quantity can be directly measured by scan- 
ning tunneling spectroscopy probes. Evaluating the sum 
over m in (|11|) . we obtain LDOS map shown in Fig[U 

Several quasi-bound states are seen in LDOS maps 
(FigH|) as local resonances at e < 0. The values of (3 
were chosen to illustrate that the width T of each reso- 
nance scales with e, while its spatial extent scales as 1/e, 
in agreement with our quasiclassical analysis and Eq. ([7]) . 

A distinct advantage of local probes, as opposed to 
transport, is that the supercritical impurities do not need 
to be a majority. In fact, it suffices to locate just one non- 
univalent impurity and perform STM imaging in vicin- 
ity. Alternatively, one can identify groups of two or three 
univalent impurities that together will act as one super- 
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FIG. 4: Spatial map of the local density of states (LDOS) 
near a charge impurity, Eq.fTTJ. The signatures of the quasi- 
bound states are resonances appearing at/3<— |atr~0 
and e < (marked by arrows). Note the localization length 
that scales inversely with e, and the linewidth proportional 
to e, as predicted by Eq.©. The intensity of the resonances 
is well in excess of the asymptotic value v(s) cc |e| at large r. 
Periodic modulation at e > with maxima at kr ~ nn is the 
standing wave oscillation [ll| (fc = e/hvF, £0 = W~ 3 hvp/ro). 



critical impurity, or even deliberately create such a group 
by inducing local charge by voltage applied to STM tip. 

The resonances © also give rise to anomalously strong 
Kondo-like effects. The striking property of the linewidth 
(O, namely, its proportionality to the energy (see Fig0]), 
indicates that the dwell time diverges at e — » 0. This 
divergence compensates the suppression of Kondo tem- 
perature by the Dirac |e| density of states. Standard 
estimates [31| for the Anderson model with the localized 
spin state associated with one of our resonant levels yield 
the Kondo temperature exponent that exhibits no sup- 
pression at £p — * 0. This is in contrast to the Kondo 



problem with extraneous spin impurities 32l. |33|. 

In summary, although massless particles are incapable 
of forming discrete states, an infinite family of quasi- 
Rydberg states can appear in a Coulomb potential of su- 
percritical strength. These quasi-bound states manifest 
themselves in a variety of physical properties, in particu- 
lar in resonant scattering and local resonances, providing 
a striking signature of the atomic collapse regime that can 
be modeled using charge impurities in graphene. 
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